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Abstract 

Here, we propose a conceptual approach for design of an ultracompact nanoscale passive optical 
circulator based on the excitation of plasmonic resonances. We study a three-port Y-junction with 
a deep subwavelength plasmonic nanorod structure integrated into its core. We show theoreti- 
cally that such a structure immersed in the magneto-optical media may function as magnetically 
tunable scatterer tilting and rotating its near-field distribution and corresponding radiation. We 
demonstrate, using numerical simulations, that such a rotation of the near-field radiation yields a 
break in the symmetry of the coupling between the junction arms and the structure in such a way 
that the signal launched from any of the three ports is mostly transmitted into the next port in 
the circular order, while the other port is essentially isolated, thus providing the functionality of 
an optical circulator with subwavelength dimensions. 
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Advances in nanophotonics and plasmonic optics in recent years have provided miniatur- 
ization of optical devices and components suitable for on-chip integration l|-|3|. In particular, 
some of the basic passive and active nanoscale optical integrated circuit elements including 
interconnects, beam splitters, switches and other components have been demonstrated re- 
cently {41. At the same time, fully functional signal processing in nanophotonics requires 
design of nonreciprocal elements, such as isolators or circulators (analogous, for example, 
with advances in microwave technology js, 6|). 

Such nonreciprocal characteristics are achieved with materials exhibiting time-reversal 
symmetry breaking, such as magneto-optic materials (MO). (Other possibilities include uti- 
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However, 



lizing nonlinearities [7H9| or materials having time dependent properties 
at optical frequencies the magneto-optical response is usually weak, which results in a bulky 
and inefficient design of basic nonreciprocal optical components, not usually compatible 
with on-chip integration. For example, "conventional" three port optical circulators - de- 



vices that guide light from any of the input ports into the next port in a circular order 12 1 

- are typically more than several hundred microns large (i.e. ^ A ^ l/im) jisl . 

Progress in silicon photonics and magneto-photonic crystals promises a substantial minia- 
turization of key nonreciprocal optical devices. For instance, an on-chip silicon-based optical 
isolation jisl, one-way edge state propagation 16, 17|, and cornpact photonic crystal circu- 
lators comparable with the operating wavelength (i.e. > A) 18l-l2l| have been suggested and 
demonstrated recently. However, further size reduction of conventional dielectric devices 
down to the subwavelength scale (< A) encounters fundamental challenges due to the light 
diffraction limit. 

This is where highly confined surface plasmons - electromagnetic fields at metal surfaces 

- can be utilized in order to extend the technology beyond the diffraction limit js], I22I . 
Furthermore, strong light localization in plasmonic structures can lead to a considerable 
increase of the nonreciprocal response. The latter was successfully implemented resulting in 
a significant enhancement of Kerr and Faraday effects 23l-l27l|. 

In our other work 28|] we have shown that plasmonic nanorod structures possessing a 
rotational symmetry, when immersed in magneto-optical media, exhibit a strong azimuthal 
symmetry breaking. In particular, it has been suggested that rotational symmetry of the 
structure gives rise to an eigenmode degeneracy at some of the plasmonic resonances. (It is 
worth mentioning here that the study of axially symmetric plasmonic nanoclusters 29l. |30| 
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and oligomers 3l| with high quahty resonances has emerged recently as a rapidly develop- 
ing and exciting field.) These eigenmodes may be represented as counter rotating eigen- 
states forming an azimuthally standing wave pattern in the unbiased regime. In presence 
of magneto-optical activity a strong mode interaction occurs yielding a frequency splitting 
of these counter-rotating states 19|] which can be described by a simple empirical formula: 
u± = Wo(l iCa); where ± correspond to two counter-rotating states, Uq is the resonant fre- 
quency at which the mode degeneracy is observed in the unbiased case, a is the strength of 
magneto-optical activity, and C is the intermodal coupling constant depending on the mode 
configuration. Note that an analogous symmetry breaking mechanism between the counter- 



rotating states was employed to design a compact photonic crystal circulator in [18|, ll9| . Our 
recent study of the light scattering from such nonreciprocal nanostrctures revealed that the 
symmetry breaking between the counter rotating states leads to a strong near-field power 



fiow circulation resulting in the tilt and rotation of the scattered radiation pattern 



In this Letter, based on our results from 
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we propose a design for an ultracompact 



nanoscale plasmonic circulator, exploiting the above-mentioned mechanism of rotation of 
scattering from subwavelength plasmonic nanorod structure in the presence of MO materials. 
We study a three-port Y-junction with a plasmonic structure embedded in its core. We show 
that the structure acts as a resonant scatterer which enhances the junction's transmission 
characteristics. We reveal that introducing MO activity into the system yields rotation of 
the near-field scattering by the subwavelength structure, and numerically demonstrate that 
such rotation breaks the symmetry between the power transmitted through the junction's 
output arms. We show that with a relatively reasonable values of MO activity the signal 
launched into any of the ports is mostly guided into the next port (in the circular order), 
while the other output port is essentially isolated. Finally, we analyze infiuence of losses 
on the system operation and show that circulator performance remains high even in lossy 
regime. 

Figure m shows the geometry of our proposed three-port junction circulator. For the sake 
of simplicity in the analysis and numerical simulations, and without loss of generality, we 
assume our system to be a two-dimensional (2D) structure. The circulator consists of three 
single-mode dielectric waveguides of width a coupled with a plasmonic structure consisting of 
metallic nanorods with deeply subwavelength diameters immersed in MO medium as shown 
in Fig. [U We assume that the MO media is magnetized along the z-direction, and therefore, 
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FIG. 1: Schematic of our subwavelength three-port junction circulator. The detailed description 
of the structure and the parameters is given in the text. 



its relative permittivity function is described by an antisymmetric tensor: 

icy Syiio 

\ £j_ y 

where Emo and e± are diagonal components, a is the off-diagonal component responsible for 
the "strength" of MO (nonreciprocal) activity of the medium. The parameter a is usually 
very small at optical frequencies (of the order of 10~^ or less). For example, bismuth iron 



garnet (BIG) exhibits Emo = 6.25, a = 0.06 and small optical losses [32|. In order to better 
understand the system dynamics of our proposed structure, in our theoretical studies we 
assume a ranging from zero to 0.15. 
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The plasmonic structure consists of a metallic nanorod with radius R = 50nm surrounded 
by three other metallic nanorods with radii Rs = lOnm as shown schematically in Fig. [TJ 
The distance between the centers of thinner and thicker nanorods is 80nm. Here we assume 
the metal relative permittivity to be fixed at = —10, and, if not stated otherwise, the 
system is assumed to be lossless. Note that, for the sake of simplicity in our analysis and 
optimization we keep the material permittivities independent of frequency. This does not 
introduce any limitation in our concept, since the dispersion of the material parameters can 
be easily included in the original design of our system. 

This 4-nanorod structure is placed in a uniform magnetized BIG host. At A ^ 1.4/im and 
the parameters chosen above, this system exhibits a sharp plasmonic resonance, a strong 
near-field power flow circulation, and rotation of the near-field scatte ring at this resonance 
(not shown here, but in full accordance with our results shown in Ref. 28|). Integrating this 
structure into the core of the three-port junction, as shown in Fig.[Tl we expect that the near- 
field circulation cause the asymmetric and nonreciprocal power redistribution into the output 
arms of the junction. In order to demonstrate our predictions we numerically study the 
circulator operation with a finite-element-based method in the COMSOL Myltiphysics"'"^. 

We begin our analysis with a magnetically inactive system, i.e. unmagnetized case when 
a = 0. In this case no asymmetry or nonreciprocity is present and the system behaves 
as a symmetrical beam splitter, i.e. signal launched into one of the ports is distributed 
symmetrically and equally between the two output ports of the junction. We search for 
the optimum regime of beam splitter operation in order to maximize the output power 
and reduce the radiation losses. For this reason we scan over the wavelength range near 
the structure resonance, and by varying the junction parameters we trace for the maximum 
power transmission. Following this procedure we find that the optimum width of waveguides 
is a = 350nm and their permittivity is 5 = 3. For matching purposes, we also introduce 
a 0.9/im-long transition waveguides with permittivity 3.5 (see Fig. [1]) in order to reduce 
backscattering and increase the transmission. We note that by placing the structure into 
the core of the junction its resonance is shifted to 1.43/im due to a change in boundary 
conditions. 

In Fig. |2]^a) we plot the simulation results for the z-oriented magnetic field distribution 
at the resonance for a signal launched into port 1. Clearly a symmetric beam splitting 
into the other two output arms is observed. In order to trace the operation efficiency we 
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FIG. 2: Simulation results for the distribution of the z-oriented magnetic field {H^) in the circulator 
at the free-space wavelength A = 1.43/im for a) unbiased (a = 0), and b),c) biased {a = 0.1) 
systems, respectively. Arrows show the direction of the power flow. Numbers 1 — 3 label the ports 
of the circulator. 

study the normalized power transmission 6*21(31) = Pouf / PL with frequency variation, see 
Fig.|4](a). We observe a strong transmission resonance at free-space wavelength A ~ 1.43/im, 
which corresponds to the plasmonic resonance of the entire structure. At the resonance, the 
power output at each of the ports 2 and 3 (5*21 = S31) reaches up-to 33% and then steeply 
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decreases with frequency detuning. The strong dependence and correlation between the 
power output and the structure excitation reveals a crucial role of structure resonance in 
the junction operation. The structure being an efficient plasmonic scatterer when excited 
redistributes the input power between the junction arms, i.e. "directionally" scatters the 
input signal and thus increasing the out-coupling efficiency. Off resonance the scatterer is 
not excited as much, and therefore the system behaves as a conventional photonic Y-j unction 
splitter, which is rather inefficient due to a large angle between the junction arms so that 
most of the input power is radiated away. We note that, despite of significant increase in 
transmission with the use of the nanostructure at its resonance, the radiation losses remain 
quite high; however we believe that further geometry optimization, including the redesign 
of the plasmonic structure itself, will increase the performance. 




FIG. 3: Snapshot in time of the near-field distribution of the magnetic field (Hz) excited at the 
resonance of the plasmonic structure (A = 1.43/xm) for a) unbiased (a = 0) and b) biased (a = 0.1) 
systems, respectively. The arrows show schematically the coupling between input port and the two 
output ports. 

As the next step, we introduce the MO activity (a ^ 0) and study its influence on the 
circulator operation. Fig. EJ^b) shows the numerical results for the z-oriented magnetic fleld 
distribution for the signal launched into the port 1 at the resonance frequency (A = 1.43/im) 
when a = 0.1. In this case we observe a substantial breaking of the symmetry between the 
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output at ports 2 and 3. The power is predominantly guided into the arm 2, while the other 
arm carries much less power. To reveal the mechanism of such strong symmetry breaking 
we study the near-field distributions for biased and unbiased systems, see Fig. [31^a,b). The 
near-field pattern in the unbiased case, Fig. ^a), is symmetric and resembles that of a 
"point" magnetic quadrupole. Such near-field distribution ensures a strong and symmetric 
coupling between the junction arms and the structure. When MO activity is introduced 



the nearfield radiation is tilted and rotated 
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This rotation of the scattering pattern 



leads to the change of coupling between the structure and the junction arms, so that input 
(port 1) and only one of the output (port 2) waveguides are coupled efficiently with the 
structure. Furthermore, we notice from Fig. El^b) that a predominantly antisymmetric field 
profile is induced at the entrance of one of the output arms (port 3 in this case), which is 
why the single-mode waveguide is not excited, causing strong isolation of the corresponding 
port, as shown on Fig. Mjo)- We note that the suggested mechanism of circulator operation 
based on the scattering and the rotation of the scattering differs conceptually from the well 
established interference principle employed in the conventional circulator designs j^. 
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FIG. 4: a) Variation of the circulator power transmission into the port 2 and isolation of the port 
3, as a function of wavelength for biased (a = 0.1) and unbiased (a = 0) cases, b) Variation of the 
circulator power transmission as a function of MO activity parameter at the free space wavelength 
(A = 1.43/im). Here metal permittivity is assumed to be frequency indedpendent and fixed at 

Em = -10. 



Next we trace the circulator performance as a function of wavelength when a = 0.1, see 
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Fig. m^a). We observe again a resonant behavior which clearly manifests the role of the 
structure in the circulator operation. At resonance we find that the power transmitted to 
port 2 comprises up-to 63% of the input power, whereas less than 3% is transmitted from 
port 1 to port 3 demonstrating a good isolation of the corresponding port. We note that 
63% of power output observed in the biased case is in fact about 96% of total power output 
in the unbiased case (33% + 33%), meaning that the radiation losses in the biased regime 
stay almost the same as in the unbiased one. 

In Fig. m^b) we plot the transmission at ports 2 and 3 as a function of MO activity 
parameter a at A = 1.43/im. We observe a growth of power output for the port 2 with a 
corresponding increase of isolation for the port 3. The transmission through port 2 peaks 
with maximum of 64% when a = 0.11 and begins to decline with a subsequent increase of 
a. For small values of a we observe linear dependence of transmission and isolation on a. 
which correlates well with the empirical formula describing the frequency splitting between 
the counter rotating states. The deviation from the linear dependence with growth of a 
can be explained by slight shift of resonance frequency and the change of coupling efficiency 
between the output waveguides and the structure due an "over rotation" of the near-field 
pattern. We note that at a = 0.06 corresponding to the realistic levels achieved in BIG we 
observe about 54% of transmission from port 1 to port 2 and less than 10% from port 1 to 
port 3. 

In order to underline the nonreciprocal nature of the proposed plasmonic circulator we 
launch signal from port 2. In this case we clearly see that light is emitted from port 3, 
whereas port 1 is isolated, as expected. Hence, the device when biased transmits the signal 
from any port to the next port (in circular order). Changing the direction of magnetization 
will swap the output and isolated ports, as also expected. Note also that the design principle 
we utilized based on the scattering from the structure provides the potential for the design 
of a circulator with higher number of ports. 

Finally, we study the infiuence of losses on the performance of the proposed circulator. 
For this purpose we consider the metal to be lossy, described by a complex permittivity 
Em = —10 — O.OH. In Fig. [5] we plot the transmission as a function of frequency for both 
output and isolated ports when a = 0.1. Clearly the symmetry breaking in transmitted 
power between the output ports is still observed. Again we achieve resonant transmission 
though one of the ports only, while the other output port is isolated. However, as expected 
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FIG. 5: The effect of metal losses on the circulator performance. Power transmission as a function 
of wavelength for a = 0.1 is plotted for the output and isolated ports. Here metal permittivity is 
assumed to be frequency indedpendent and fixed at = — 10 — O.Oli. 



the losses result in decrease in the total power transmission, and a decrease in the isolation 
ratio, i.e. the ratio of transmitted powers through both output ports. The ability to keep 
the circulator performance with losses allows the possibility of implementing ferromagnetic 
metals possessing stronger MO activity (such as, for example, Co, Ni, Fe). Although these 
materials exhibit quite high electromagnetic damping, their combination with plasmonic 
nobel metals and proper structure engineering might give the possibility to control the 



effects of the losses and the MO activity 2J, |25 |. 



In conclusion, we have proposed a design for a subwavelength three-port junction circu- 
lator with an integrated plasmonic structure. We have demonstrated numerically that the 
plasmonic resonance and the rotation of corresponding near-field scattering in presence of 
MO activity leads to the breaking of symmetry in the power output between the circulator 
arms. Our analysis reveals up-to 63% of power transmitted through one of the output ports, 
with almost complete isolation for the other output port of the circulator. Finally, we have 
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studied the influence of losses and found that despite the decrease of net output in lossy 
systems, the symmetry breaking between the output ports remains relatively high. This 
concept may find useful application in integrated nanophotonics, optical metatronics, and 
signal handling in on-chip integrated circuits. 

I. ACKNOWLEDGEMENTS 

Authors thank A. Mahmoud for useful discussion. This work is supported in part by the 
US Air Force Office of Scientific Research (AFOSR) grant number FA9550-10-1-0408. 



[1] Murphy, E.J., Ed. Integrated Optical Circuits and Components: Design and Applications; 

Marcel Deckker: New York, 1999. 
[2] Joannopoulos, J.D.; Villeneuve, P.R.; Fan, S. Nature 1997, 386, 143-149. 
[3] Bozhevolnyi, S.I., Ed. Plasmonic Nanoguides and Circuits; Pan Stanford: Singapore, 2009. 
[4] Ebbesen, T.W.; Genet, C; Bozhevolnyi, S.I. Phys. Today 2008, 61, 44-50. 
[5] Helszajn, J. Nonreciprocal microwave junctions and circulator, Wiley, 1975. 
[6] Helszajn, J. Microwave Engineering: Passive, Active and N on- Reciprocal Circuits; Mcgraw- 

Hih, 1992. 

[7] Shadrivov, I.V.; Fedotov, V.A.; Poweh, D.A.; Kivshar, Y.S.; Zheludev, N.I. New J. Phys. 
2011, 13, 033025-8. 

[8] Miroshnichenko, A.E.; Brasselet, A.E.; Kivshar, Y.S. Appl. Phys. Lett. 2010, 96, 063302-3. 
[9] Alberucci, A.; Assanto, G. Opt. Lett. 2008, 33, 1641-1643. 
[10] Yu, Z.; Fan, S. Nature Photon. 2009, 3, 91-94. 

[11] Lira, H.; Yu, Z.; Fan, S.; Lipson, M. Phys. Rev. Lett. 2012, 109, 033901-5. 

[12] Fay C.E.; Comstock, R.L. IEEE Trans. Microwave Theory Tech. 1965, 13, 15-27. 

[13] Takei, R,; Mizumoto, T. Japan. J. Appl. Phys. 2010, 49, 052203-6. 

[14] Sugimoto, N.; Shintaku, T.; Tate, A.; Terui, H.; Shimokozono, M.; Kubota, E.; Ishii, M.; 

Inoue, Y. IEEE Photon. Tech. Lett 1999, 11, 355-357. 
[15] Bi, L.; Hu, J.; Jiang, P.; Kim, D.H.; Dionne, G.F.; Kimerling, L.C.; Ross, C.A. Nature Photon. 

2011, 5, 758-762. 



11 



[16] Wang, Z.; Chong, Y.; Joannopoulos, J.D.; Soljacic, M. Nature 2009, 461, 772-775. 

[17] Liu, S.; Lu, W.; Lin, Z.; Chui, S.T. Appl. Phys. Lett. 2010, 97, 201113-3. 

[18] Wang, Z.; Fan, S. Opt. Lett. 2005, 30, 1989-1991. 

[19] Fan, S.; Wang, Z. J. Magn. Soc. Jpn. 2006, 30, 641-645. 

[20] Smigaj, W.; Romero- Vivas, J.; Gralak, B.; Magdenko, L.; Dagens, B.; Vanwolleghem, M. Opt. 

Lett. 2010, 35, 568-570. 
[21] Qiu, W.; Wang, Z.; Soljacic, M. Opt. Express 2011, 19, 22248-22257. 
[22] Maier, S.A. Plasmonics: Fundamentals and Applications; Springer- Verlag: Berlin, 2007. 
[23] Belotelov, V.I.; Akimov, LA.; Pohl, M.; Kotov, V.A.; Kasture, S.; Vengurlekar, A.S.; Gopal, 

A.V.; Yakovlev, D.R.; Zvezdin, A.K.; Bayer, M. Nature Nanotech. 2011, 6, 370-376. 
[24] Wang, L.; Clavero, C; Huba, Z.; Carroll, K.J.; Carpenter, E.E.; Gu, D.; Lukaszew, R.A. Nana 

Lett. 2011, 11, 1237-1240. 
[25] Jain, P.K.; Xiao, Y.; Walsworth, R.; Cohen, A.E. Nano Lett. 2009, 9, 1644-1650. 
[26] Tomita, S.; Kato, T.; Tsunashima, S.; Iwata, S.; Fujii, M.; Hayashi, S. Phys. Rev. Lett. 2006, 

96, 167402-4. 

[27] Sepulveda, B.; Gonzalez-Dyaz, J.B.; Garcya-Martyn, A.; Lechuga, L.M.; Armelles, G. Phys. 

Rev. Lett. 2010, 104, 147401-4. 
[28] Davoyan, A.R.; Engheta, N. \arXiv: 1212.6585 2012. 

[29] Fan, J. A.; Wu, C; Bao, K.; Bao, J.; Bardhan, R.; Halas, N.J.; Manoharan, V.N.; Nordlander, 

R; Shvets, G.; Capasso, F. Science 2010, 328, 1135-1138. 
[30] Lassiter, J.B.; Sobhani, H.; Fan, J. A.; Kundu, J.; Capasso, F.; Nordlander, P.; Halas, N.J. 

Nano Lett. 2010, 10, 3184-3189. 
[31] Hentschel, M.; Dregely, D.; Vogelgesang, R.; Giessen, H.; Liu, N. ACS Nano 2011, 5, 2042- 

2050. 

[32] Tepper, T.; Ross, C.A. J. Cryst. Growth 2003, 255, 324-331. 



12 



